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computer  run  time. 

•  The  trial  particle  size  group  activity  distribution  reduces  the  size  of  the 
close-in  contours  and  allows  a  slight  improvement  to  the  SEER/CIVIC 
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SUMMARY 


A  trial  methodology  was  developed  for  the  SEER  fallout  code  In  order  to 
conserve  total  fallout  activity  and  Improve  tiie  approximation  to  DELFIC  results.  The 
trial  methodology  was  compared  with  DELFIC  and  also  with  the  CIVIC  version  of 
SEER,  which  utilizes  an  adjusted  K-factor  to  improve  the  approximation  to  DELFIC. 
Also,  the  particle  size  group  activity  distribution  currently  used  in  SEER  was  replaced 
by  one  approximately  the  same  as  that  used  in  DELFIC.  Combinations  of  trial  and  old 
methodology  with  trial  and  old  activity  distributions  were  adl  run  for  comparison  with 
DELFIC.  The  results  of  the  comparison  show  that: 

•  The  trial  methodology  conserves  activity  while  utilizing  the  original  K- 
factor. 

•  The  trial  methodology  produces  fallout  contours  which  are  a  better 
approximation  to  DELFIC  than  are  those  produced  by  the  old  SEER  with 
the  original  K-factor. 

•  The  current  CIVIC  version  of  SEER  with  the  adjusted  K-factor  is  a  slightly 
better  approximation  to  DELFIC  than  the  trial  methodology. 

•  The  trial  methodology  requires  on  the  average,  about  20%  additional 
computer  run  time. 

•  The  new  particle  size  group  activity  distribution  reduces  the  size  of  the 
close-in  contours  and  allows  a  slight  improvement  to  the  SEER/CIVIC 
approximation  to  DELFIC  at  the  3000  r/hr  contour,  and  a  slightly  less 
accurate  approximation  at  the  100  and  300  r/hr  contours. 


It  is  possible  that  further  work  on  the  trial  methodology  could  improve  the 
approximation  to  DELFIC.  Pending  this,  however,  it  is  concluded  from  this  effort  and 
the  results  of  Reference  1  that  the  current  CIVIC  version  of  SEER  with  adjusted  K- 


factor  is  the  best  available  fast-running  approximation  to  DELFIC. 
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SECTION  1 


INTRODUCTION 

This  report  presents  the  methodology  and  results  of  an  attempt  to  modify 

and  Improve  the  computer  code  SEER  (Simplified  Estimation  of  Exposure  to  Radia- 

•-  — 

tion).  SEER  ist  part  of  the  large  nuclear  weapon  damage  assessment  code  CIVIC. 
.SEER-(Relerences  2,  ^  was  developed  in  the  early  1970s  with  the  intent  that  it 
simulate  the  fallout  exposure  rate  contours  of  DELFIC  (Defense  Land  Fallout 
Interpretive  code)  while  using  a  minimal  amount  of  computer  time. 

A  problem  that  surfaced  with  the  development  of  SEER  was  the  fact  that 
the  SEER  methodology  (see  S^tion  2)  leads  to  an  overprediction  of  the  total  activity 
deposited  on  the  ground  (see  References  4,  5).  In  order  to  compensate  for  this 
characteristic,  SAI  has  lowered  the  value  of  the  K-factor  in'  the  SEER  portion  of  CIVIC 
to  approximately  half  of  the  original  value.  The  K-factor  is^the  ajnount  of  weapon 
fallout  activity  per  unit  yield,  distributed  over  unit  area  (see  Section  4-1.3)i  Its  value 
is  not  precisely  known,  in  part  because  the  percentage  of  total  radioactive  debris  that 
constitutes  early  fallout  (as  opposed  to  that  which  is  lofted  to  high  altitudes  and 
carried  off  beyond  local  fallout  distances)  is  not  well  known  and  varies  with  height  of 
burst  and  other  conditions.  The  original  K-factor  used  in  SEER  (Reference  3,  p.  15)  is 

I 

tciken  from  Reference  6.  ,  The  value  indicated  by  Reference  7  is  about  25%  lower  than 
this. 

This  work  effort  investigated  the  possibility  of  modifying  the  SEER 
methodology  so  that  activity  is  conserved  and  a  good  approximation  to  DELFIC  results 
could  be  obtained  without  resort  to  the  adjustment  of  the  K-factor.  The  following 
sections  describe  the  current  SEER  methodology,  the  trial  methodology,  and  the 
results  of  the  comparison. 
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SECTION  2 


CURRENT  SEER  APPROACH 

As  shown  in  Figure  2-1,  SEER  considers  a  cylindrical  stabilized  cloud  and 
divides  it  into  (NLVL-1)  slabs,  the  number  depending  on  the  weapon  yield.  The  cloud  is 
further  discretized  by  dividing  the  activity  among  25  particle  groups  (subscript  IB), 
each  of  which  is  assigned  a  fraction  of  the  total  activity  as  well  as  characteristic  fall 
times  from  various  altitudes.  The  assumption  is  made  that  each  altitude  slab  has  an 
equal  share  of  the  total  activity,  and  that  the  particle  distribution  is  the  same  in  each 
altitude  slab.  SEER  calculates  the  ground  locations  and  radii  of  25  x  NLVL  discs,  each 
disc  orginating  at  one  of  the  discreet  slab  boundary  altitudes  corresponding  to  IA=  1,2, 
...NLVL.  The  situation  is  shown  in  Figure  2-2,  where  the  A,B,C,...  "hot  lines" 
correspond  to  altitude  levels  lA  =  1,2,3...,  and  the  disc  landing  points  1, 2,3,4. ..corre¬ 
spond  to  the  various  particle  size  groups  IB  =  1,2,3,...25.  SEER  computes  the  landing 
points  of  the  centers  of  each  disc.  The  activity  densities  of  the  discs  are  calculated 
using  an  area  spread  which  is  a  multiple  of  the  disc  diameter  and  a  fraction  of  the 
distance  from  the  GZ  or  from  the  previous  disc  (SEER  3.51-SEER  3.92).  Subsequent 
activity  densities  are  calculated  using  a  log  smoothing  procedure.  Following  these 
calculations  is  an  entire  set  of  calculations  which  compute  the  fallout  densities  at  a 
set  of  "key  points",  as  illustrated  in  Figure  2-3.  The  rationale  stated  for  this 
additional  step  is...  "so  that  an  appropriate  network  of  exposure-rate  values  will  be 
available  for  the  mapping  routines."  Consideration  of  contributions  from  more  than 
one  disc  to  each  key  point  is  limited  solely  to  discs  from  adjacent  altitude-level  "hot 
lines"  that  lie  within  a  disc  radius  of  the  key  point.  Moreover,  this  is  done  only  for 
discs  of  the  same  particle  class.  As  might  be  suggested  by  Figure  2-2,  this  whole 
procedure  can  lead  to  "gaps"  of  no  activity  density  between  discs.  Overestimation  of 
total  fallout  activity  occurs  when  the  activity  of  map  points  that  fall  in  these  gaps  is 
"filled  in"  through  the  SEER  process  of  interpolating  between  key  points  to  arrive  at 
map  point  activity  densities.  In  summary,  the  possible  causes  of  inaccuracies  in  SEER 
are  as  follows: 

•  Initial  use  of  an  incorrect  area  (SEER  3.79)  spread  (or  "smear")  of  activity 
from  one  disc,  based  only  on  the  diameter  of  the  disc  itself  and  its  distance 
from  the  GZ  or  the  previous  disc  (SEER  3.51-SEER  3.92). 


•  PARTICLE  CLASS  LANDING  POINTS 

o  PATTERN  EDGE  POINTS 


Note;  Paths  A,  B,  and  C  are  the  particle  class  paths  for  3  cloud  levels. 
The  example  shows  only  4  particle  classes  but  the  SEER  model 
uses  25  particle  classes.  The  L  and  R  points  are  along  the  left 
and  right  edges  of  the  fallout  pattern. 


Figure  2-2.  Example  of  landing  pattern  for  the  SEER  particle  classes 


PARTICLE  CLASS  LANDING  POINTS 
PATTERN  EDGE  POINTS 
FALLOUT  PATTERN  KEY  POINTS 


Figure  2-3.  Location  of  key  points  in  the  SEER  fallout  pattern 


•  Failure  to  consider  contributions  of  discs  from  other  particle  size  groups  to 
the  activity  at  the  given  point  (either  disc  impact  point  or  "key  point") 
(SEER  3). 

•  Failure  to  consider  contributions  of  discs  from  other  altitude  levels  beyond 
disc  radius  distance  (SEER  3). 

•  "Filling  in"  of  activity  by  interpolation  at  locations  where  the  activity  may 
not  exist  (or  be  less  than  the  interpolated  value). 


It  is  also  true  that  SEER  does  not  include  a  model  for  the  cloud  stem.  This 
would  tend  to  lower  the  amount  of  close-in  fallout,  but  SEER  compensates  (to  a 
certain  degree)  for  this  by  utilizing  a  particle  size  group  activity  distribution  which 
places  more  activity  in  the  largest  particle  size  groups  than  does  DELFIC. 


SECTION  3 


MODIFIED  APPROACH 

In  view  of  the  considerations  discussed  in  the  previous  section,  an 
alternative  approach  is  offered  which  will  eliminate  the  entire  SEER  procedure 
following  SEER  2,  which  calculates  the  positions  and  radii  of  the  25  x  NLVL  discs.  The 
procedure  which  is  substituted  is  discussed  in  the  following  sections.  In  brief,  rather 
than  establish  intensities  at  arbitrary  "key  points"  and  determine  map  point  intensities 
by  interpolation,  the  substitute  approach  "smears"  the  activity  from  each  disc  over  an 
area  corresponding  to  altitude  levels  lA  to  IA+1  and  particle  size  groups  IB-)4  to  IB+)4. 
Each  quadrangular  area  (actually  a  "rounded  off"  quadrangle)  is  then  investigated  to 
determine  which  map  grid  points  fall  within  it,  and  the  intensity  of  the  map  grid  point 
is  augmented  accordingly. 

3-1  SPREAD  OF  FALLOUT  ACTIVITY  FROM  ONE  DISC  (SUBROUTINE  SMEAR) 

Consideration  of  the  area  over  which  the  activity  from  one  disc  (say  IB,IA) 
should  be  spread  can  be  facilitated  by  referring  to  Figure  3-1.  SEER  2  calculates  the 
positions  XTAY  (IB,IA),  YTAY(IB,IA)  for  1<IB<25  and  1<IA<NLVL.  The  disc  radius 
RDTAY  (IB,  lA)  is  also  calculated.  The  disc  initial  altitude  is  at  the  bottom  of  the 
finite  altitude  slab  that  it  represents  (see  Figure  2-1).  Thus,  the  activity  from  that 
slab  is  not  confined  to  the  IB,  lA  disc,  but  rather  is  distributed  between  that  disc  and 
the  disc,  IB,  lA  +  1,  representing  the  next  higher  altitude  slab.  Moreover,  the  particle 
size  group  represented  by  the  subscript  IB  is  actually  a  range  of  particle  sizes  (e.g., 
IB=1  represents  particles  diameters  in  the  range  2241-8768  microns).  The  fall  times 
(read  in  as  data  statements  in  subroutine  BLKONE)  represent  an  average  value  over 
the  size  group.  Therefore  the  IB,  lA  disc  activity  should  also  be  spread  over  an  area 
between  the  adjacent  particle  size  groups  IB-1  and  IB  +1.  We  thus  arrive  at  an 
effective  area  "smear"  of  activity  represented  in  Figure  3-2.  The  proper  spread  is 
enclosed  in  quadrangle  PI 3,  P23,  P33,  P43,  with  the  hatched  corner  areas  excluded.  It 
might  appear  from  Figure  3-2  that  the  use  of  this  "outer  quadrangle"  would  create  an 
overlap  strip  with  artifically  high  activity  density.  In  reality,  however.  Figure  3-2 


r(iB*i,  lA) 
niB+l,  lA) 


XTAY{1B+1,  IA+1) 
YTAY(IB+1,  lA+l) 
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XTAY( IB-1,  lA) 
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Figure  3-1.  Represent.'Ytlve  SEER  results  following  subroutine  SEER  2. 
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represents  only  the  situation  for  the  very  lightest  of  the  particle  groups.  For  most  of 
the  groups,  there  is  considerable  overlap  of  these  "quadrangles"  because  the  distances 
Si,  S2,  S3,  Sk  are  less  (in  some  cases,  much  less)  than  the  disc  diameters.  Even  at  the 
outer  regions,  where  the  lightest  particles  play  a  part,  the  overlap  does  not  appear  to 
present  a  problem.  Indeed,  the  modified  SEER  provides  a  pattern  which  correlates 
well  with  that  of  DELFIC. 

3-1.1  Finding  the  Pbints  Pi(XiYi),  P2(X2,Y2),  P3(X3tY3),  and  P4(X4,Y4) 

These  points  are  approximated  as  midpoints  between  adjacent  particle 
group  wafer  locations.  Thus: 

Xl=  (XTAY(IB-1,IA)  +  XTAY  (IB,IA))/2. 

X2  =  (XTAY(IB,IA)  +  XTAY(IB+l,IA))/2. 

...  etc. 

However,  if  IB=1,  we  approximate  as  follows: 

XI  =  XTAY  (IB,  lA)-  DELTX 

Where  DELTX  =  X2  -  XTAY(IB,IA) 

...X2  having  been  calculated  first. 

If  IB=23,  we  approximate  as  follows: 

X2  =  XTAY(IB,IA)  +  DELTX 
where  DELTX  =  (XTAY(lB,IA)-XTAY(IB-l,IA))/2 

Analagous  approximations  are  done  for  points  P3  and  P^^.  Note  that  we  do 
not  have  to  make  interpolation  approximations  for  altitude  increments.  SEER  has 
already  calculated  the  location  of  the  discs  for  1A=NLVL,  which  represents  the  upper 
limit  for  the  top  slab  NSLABS=NLVL-1. 
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3-1.2  Calculating  the  Average  Wafer  Radius,  R 

Inspection  of  SEER  output  shows  that  over  most  of  the  particle  size  groups, 
the  value  of  RDTAY  (IB,IA)  does  not  change  drastically  from  one  disc  to  the  next. 
Thus,  in  the  interest  of  conserving  computer  time,  it  is  deemed  sufficiently  accurate 

to  represent  the  disc  radii  corresponding  to  each  of  the  points  PpP2’  ^3’  ^4'  ^ 

single  disc  size  with  a  radius  which  is  the  average  over  the  six  discs  shown  in  Figure  4- 
1,  with  discs  IB,IA  and  IB,  IA+1  each  weighted  twice. 

3-1.3  Calculate  the  Cross-Ehxxhicts  of  Quadrangle  PiP2P3P% 

In  order  to  locate  map  points  in  or  out  of  the  ”rounded-off”  quadrangle,  it 
will  be  necessary  to  compare  cross-product  signs  (see  Section  3-3.2  below).  It  is 
therefore  necessary  to  obtain  the  cross  products  of  each  side  and  its  adjacent  side. 
Since  it  is  only  the  signs  of  the  cross-products,  not  their  value,  that  are  important,  it 
is  sufficient  to  obtain  the  cross-products  of  the  "inner”  quandrangle  (PI,  P2,  P3,  P4). 
As  an  example  to  illustrate  the  method  for  obtaining  the  cross-products,  we  will 
consider  ^  x  For  a  two-dimensional  cross-product  (or  vector-product)  in  the  x-y 
plane,  the  cross  product  is  obtained  from: 

^  X  ^  =  (ajb2-a2bj)k 

where  k  is  the  unit  vector  in  the  2-direction  (into  or  out-of  the  page),  and  the  vectors 
are  represented  by; 

+bji 

S2  =  a2i.+  ^21- 

where  i_and  j_are  the  unit  vectors  in  the  X  and  Y  directions,  respectively. 

Now,  aj  =  X2  -XI  =  DXl  ( in  subroutine  SMEAR) 

bj=  Y2  -  Y1  =  DYl 


so  M  X  S2  =  DX 1  *  DY2  -  DX2  *  DY 1  =  CP1X2  ( in  subroutine  SMEAR) 

likewise  CP2X3  =  DX2  *  DY3  -  DX3  *  DY2 

CP3X4  =  DX3  *  DY4  -  DX4  *  DY3 
CP  4X1  =  DX4  *  DYl  -  DXl  *  DY4 

The  cross-products  are  checked  to  see  if  they  are  .GT.O  or  .LT.O  and  if  they  all  agree 
in  sign.  If  they  do  not  all  agree  in  sign,  then  an  error  message  is  printed. 

3-1.4  Locating  the  Tangent  Points  Pll,  P12,  etc. 

In  order  to  construct  the  "outer"  quadrangle  of  sides  S^'  $2'  S^'  S^',  the 
tangent  points  Pll,  P12,  P21,  P22,  P31,  P32,  P41,  P42  are  first  located  (see  Figure 
3-2).  In  order  to  locate  these  points,  it  is  first  necessary  to  determine  if  the  inner 
quadrangle,  defined  by  points  PI  P2  P3  P4  in  that  order  ,  progresses  in  a  clockwise  or 
counterclockwise  direction  (the  reason  for  this  is  discussed  in  the  next  section).  If  the 
progression  from  PI  to  P2  to  P3  to  P4  is  clockwise,  then  the  cross-products  of  the 
sides  of  the  inner  quadrangle  will  be  negative.  If  counterclockwise,  then  the  cross- 
products  will  be  positive.  The  cross  products  discussed  here  are  as  follows:  ^  x  S2,  ^ 
X  ^  S3  X  ^  and  S4  x  SI,  and  are  obtained  as  described  in  Section  3-1.3. 

As  an  example  for  this  procedure,  we  consider  Figures  3-3  and  3-4.  In 
Figure  3-3,  we  have  assumed  for  example  that  the  inner  quadrangle  progression  PI  P2 
P3  P4  is  clockwise.  We  know  the  locations  of  PI  and  P2  and  it  is  apparent  that  the 
polygon  PI  P12  P2I  P2  is  a  rectangle.  Therefore,  the  slope  of  line  SRI  is  the  negative 
inverse  of  the  slope  (SMI)  of  line  SI.  We  can  therefore  write  2  equations: 

YRl  /  XRl  =  -1./SM1 

YRl^  +  XRl^  =  SRl^  Where  SRI  =  R  in  magnitude 

Solving  these  simultaneously: 

X12  =  X1-FX*R 
Y12  =  Y1  +  FY  »  R 

FX  =  l./SQRT(l.  +  1./SM1  **  2) 

FY  =  I./SQRT(1.  +  SMI  **  2) 
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This  particular  example  has  been  case  3  of  the  clockwise  case.  To  understand  its 
context  in  the  overall  calculation,  and  the  other  equations  involved,  consider  Figure 
3-4  and  Table  3-1.  Figure  3-4  shows  the  various  possible  orientations  of  the  lines  SI 
(hence  SI'),  S2  (hence  S2'),  S3  (hence  S3'),  and  S4  (hence  S4')  along  with  the 
corresponding  "case  numbers".  Table  3-1  shows  the  resulting  equations  for  the  X  and 
Y  locations  of  the  points  PI  2  and  P21  for  each  case  in  both  the  clockwise  and 
counterclockwise  conditions.  We  can  see  from  Figures  3-3  and  3-4  and  Table  3-1  the 
necessity  of  knowing  the  clockwise/counterclockwise  condition  (hence  the  cross- 
product  signs).  For  the  CCW  case,  the  tangent  points  all  lie  to  the  right  of  the  points 
PI,  P2,  P3,  P4.  For  the  CW  case,  they  all  lie  to  the  left. 

3-1.5  Finding  the  Outer  Quadrangle  Vertices  P13,  P23,  P33,  P43 

The  procedure  for  finding  these  points  is  to  find  the  general  equations  for 
the  lines  SI',  S2',  S3',  and  S4'  and  solve  each  intersecting  pair  simultaneously  to  locate 
the  point  of  intersection.  The  general  form  of  the  equation  of  a  line  is: 

AX  +  BY  >  C  =  0 

Let  us  look  at  the  particular  point  PI 3,  and  the  equations  that  will  lead  to  its  location. 
The  slopes  of  the  lines  SI'  and  S4'  are  equal  to  those  of  SI  and  S4;  i.e.,  SMI  and  SM4. 
We  can  therefore  write; 

Y-Y12  =  SM1*(X-X12) 

Y-Y42  =  SM4»(X-X42) 

rearranging  to  get  the  general  form,  we  get; 

Y-SMUX  +  (SMUX12-Y12)  =  0 
Y-SM4*X  +  (SM4*X42-Y42)  =  0 


solving  for  X: 


X13(IB,IA)  =  (X42»SM4-X12*SMUY12-Y42)/(SM4-SM1) 


Table  3-1.  Equations  for  the  tangent  points  P12,  P21  for  all  orientations  of  SI 
(equations  for  other  tangent  points  are  analogous  for  other  lines  S2,  S3,  S^). 


Case 

Number 

1 


2 


3 


4 


5 


6 


7 


S 


Slope 

SMI 

0 


0 


.GT.O 


.GT.O 


.LT.O 


.LT.O 


oo 


X2-X1 


=  DXl 

CW  Equations 

CCW  Equations 

.GT.O 

X12=X1 

Y12=Y1+R 

X2UX2 

Y2UY2+R 

X12=X1 

Y12=Y1-R 

X2UX2 

Y2UY2-R 

.LT.O 

X12=X1 

Y12=Y1-R 

X21=X2 

Y2UY2-R 

X12=X1 

Y12=Y1+R 

X21=X2 

Y21=Y2+R 

.GT.O 

X12=X1-FX*R 

Y12=Y1+FY*R 

X21=X2-FX*R 

Y2UY2+FX»R 

X12=X1+FX»R 

Y12=Y1-FY*R 

X21=X2+FX*R 

Y21=Y2-FY*R 

.LT.O 

X12=X1+FX*R 

Y12=Y1-FY*R 

X2l=X2+FX*R 

Y2l=Y2-FY*R 

X12=X1-FX»R 

Y12=YUFY*R 

X21=X2-FX*R 

Y21=Y2+FY»R 

.GT.O 

X12=X1+FX*R 

Y12=Yl+FY*R 

X21=X2+FX*R 

Y21=Y2+FY»R 

X12=X1-FX*R 

Y12=Y1-FY*R 

X21=X2-FX*R 

Y21=Y2-FY*R 

.LT.O 

X12=X1-FX*R 

Y12=Y1-FY»R 

X21=X2-FX*R 

Y2i=Y2-FY*R 

X12=X1+FX*R 

Y12=YUFY*R 

X21=X2+FX*R 

Y21=Y2+FY*R 

DYl.GT.O 

X12=X1-R 

Y12=Y1 

X21=X2-R 

Y21=Y2 

X12=X1+R 

Y12=Y1 

X2UX2+R 

Y2UY2 

DYl.L.T.O 

X12=XUR 

Y12=YI 

X21=X2+R 

Y21=Y2 

X12=X1-R 

Y12=Y1 

X21=X2-R 

Y21=Y2 

OO 


If  we  solve  for  Y,  we  get: 

Y13(IB,IA)  =  (Y12*SM4-Y42*SM1+(X42-X12)*SM<**SMI)/(SM4-SM1) 

The  equations  for  P23,  P33,  and  P43  are  obtained  through  the  same  method. 

3-1.6  Finding  the  Total  Area  of  Activity  Spread 

The  total  area  over  which  the  activity  from  disc  (IB,  lA)  is  spread  can  be 
calculated  by  considering  it  to  consist  of  three  parts  (see  Figure  3-2); 

•  The  area  of  the  "inner"  quadrangle  PI,  P2,  P3,  P4  (AREA  0  in  Subroutine 
SMEAR). 

This  area  is  obtained  from  the  expression  for  the  area  of  a  polygon  from 
ansdytic  geometry,  and  is: 

AREA  0  (IB,IA)  =  0.5*(X1*Y2  +  X2*Y3  +  X3*Y4  +  X4*Y1 
-X2*Y1  -  X3*Y2  -  X4*Y3  -  XUY4) 

where  all  the  points  in  the  parentheses  are  subscripted  (IB,  lA). 

•  The  areas  of  the  4  quadrangles  (PI,  P2,  P21,  P12),  (P2,  P3,  P31,  P22),  (P3, 
P4,  P41,  P32),  and  (P4,  PI,  Pll,  P42).  The  total  of  these  4  areas  is  equal  to 
R*(51+S2+S3+S4). 

•  The  remaining  4  areas  from  the  corner  discs  (e.g.,  P21-P2-P22^P21).  It 
can  be  shown  that,  regardless  of  the  relative  orientations  of  the  points  PI, 
P2,  P3,  P4,  these  4  areas  total  up  to  FIR^. 

3-1.7  A  Word  About  the  Quadrangle  Comers 

The  SMEAR  theory  considers  a  spread  of  extremely  thin  discs  across  the  area 
described  above,  each  representing  an  extremely  small  increment  of  altitude  and  particle 
size.  The  extremities  of  the  ctrea  spread  will  be  the  corner  discs,  and  do  not  include  the 
corners  (the  shaded  portion  in  Figure  3-5).  However,  the  method  of  determining  if  map 
points  are  within  the  area  (Subroutine  MAPDEN  discussed  in  Section  4.3)  is  such  that  map 
points  that  lie  within  the  shaded  area  will  be  counted  unless  otherwise  excluded.  Such  an 
exclusion  routine  is  incorporated  in  Subroutine  SMEAR. 


3-1.8  Computing  the  Fallout  Activity  Density 


The  total  fallout  activity  input  to  the  modified  SEER  can  be  in  terms  of 
curies/KT  at  1  hour,  or  R/hr  at  1  hour  3  ft  above  an  infinite  plane,  or  some  other  value 
(represented  by  SEED)  input  by  the  user.  The  activity  in  each  altitude  slab  of  the 
stabilized  cloud  is: 

ACTLYR  =  SEED  *  YLD  *  FFR  *  PCT  /  FLOAT  (NSLABS) 


where 

YLD  =  yield  in  KT,  and  PCT  and  FFR  are  correction  factors  already  in  SEER. 
The  activity  density  in  the  quadrangle  is  thus 

DENS  (IB,IA)  =  ACTLYR  »  PERCNT  {IB)/AREA  {1B,IA) 

where  PERCNT  (IB)  is  the  percent  of  total  activity  ascribed  to  the  particular  IB^^ 
particle  size  group. 

3-2  SETTING  UP  THE  FALLOUT  MAP  (SUBROUTINE  GRID) 

The  purpose  in  establishing  the  map  grid  at  this  point  is  that  the  methodology 
used  in  subroutine  MAPDEN  (description  below)  requires  that  the  coordinates  of  the  map 
points  be  known.  This  is  in  contrast  with  the  current  SEER  method,  which  computes 
activity  densities  at  "key  points"  before  the  grid  is  established,  and  later  computes  the 
activity  densities  at  the  map  points  by  interpolation  between  the  key  points.  The  map 
points  established  in  GRID  are  all  located  in  terms  of  X  and  Y  distances  in  meters  from 
GZ.  The  basic  approach  to  map  construction  is  the  same  as  SEER.  That  is,  the  map  is 
printed  on  a  line  printer  with  X-values  increasing  from  left-to-right  and  y-values  printed 
in  decreasing  order  from  top-to-bottom.  The  Y-length  may  be  as  long  as  the  user  wishes 
(AMGY),  and  the  X-width  may  be  divided  into  strips  so  that  the  printer  can  handle  a  wide 
map  that  still  retains  some  reasonable  resolution.  The  user  puts  in  the  number  of  X  points 
(NCOL)  to  be  used  for  each  strip,  and  also  puts  in  an  estimate  of  the  total  number  of  X 
points  (AMGX)  required,  usually  some  multiple  of  NCOL.  An  alternate  method  is 
available  by  setting  an  index  (IMAP)  to  1.  This  allows  the  user  to  input  the  extent  of  the 
map  in  both  the  X  and  Y  directions,  and  also  the  distance  (DEL)  from  one  grid  point  to  the 
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next  (DEL  is  the  same  in  the  X  and  Y  directions).  This  feature  duplicates  a  feature  that  is 
available  in  DELFIC,and  is  particularly  useful  if  one  has  an  idea  of  the  extent  of  activity 
levels  of  interest  (say  from  a  prior  execution  of  SEER)  that  would  thus  allow  the  map 
routine  to  avoid  calculations  over  extended  areas  of  low  activity.  The  use  of  this  feature 
allows  the  program  to  skip  over  the  calculations  described  in  the  next  two  subsections. 

3-2.1  Finding  the  Max/Min  Limits  of  the  Fallout  Distribution 

This  is  achieved  by  finding  the  maximum  and  minimum  values  of  X  and  Y  for 
the  quadrangle  vertices  for  all  lA  levels  and  IB  particle  size  groups.  The  extreme  values 
found  for  each  are  called  XMIN,  XMAX,  YMIN,  YMAX. 

3-2.2  Setting  up  the  Grid 

A  grid  is  set  up  as  follows:  First  the  extents  of  fallout  in  the  X  and  Y 
directions  are  determined 

DELX  =  XMAX  -  XMIN 
DELY  =  YMAX  -  YMIN 

Next,  the  grid  increment  is  determined.  A  square  grid  increment  is  used}  that 
is,  the  increment  (DEL  )  from  one  X-point  to  the  next  is  the  same  as  that  from  one  Y-point 
to  the  next.  The  following  quantities  are  computed; 

FI  =  DELX/AMGX 
F2  =  DELY/AMGY 

and  the  quantity  F3  =  AMAX1(F1,F2)  is  taken  as  the  first  approximation  to  DEL.  F3  is 
then  increased  slightly  up  to  some  multiple  of  10,  100,  1000,  etc  (this  makes  the  map 
easier  to  read  than  the  current  SEER).  The  map  starting  points  are  established  at 
some  X  value  just  less  than  XMIN  that  is  divisible  by  DEL,  and  some  Y-value  just 
greater  than  YMAX  which  is  divisible  by  DEL.  The  number  of  X-points  in  the  map  is 
calculated  as 

NXPTS  =  IFIX  (DELX/DEL)  +  3 


Vi- 

where  3  has  been  added  to  correct  for  truncation  and  to  provide  a  margin  on  either 
side  of  the  map.  The  number  of  Y  points  is  detennined  in  a  similar  manner.  The 
number  of  map  strips  required  is  determined  by  comparing  NXPTS  with  NCOL. 

3-3  CALCULATING  THE  FALLOUT  ACTIVITY  DENSITY  AT  EACH  MAP  POINT 
(SUB-ROUTINE  MAPDEN) 

MAPDEN  loops  through  each  disc  AREA  (1B,IA)  (1A=1,  NSLABS),  (IB=1,  25) 
to  determine  which  map  points  fall  within  the  area.  This  could  require  excessive 
computer  time  if  each  IA,1B  iteration  had  to  consider  all  of  the  map  points.  In  order 
to  save  time,  the  number  of  map  points  can  be  limited  to  those  that  are  in  the  near 
vicinity  of  the  quadrangle  because  the  map  grid  point  locations  have  already  been 
established  in  Subroutine  GRID. 

3-3.1  Find  the  Limit  "Box"  for  Map  Points 

The  limits  for  map  points  to  be  considered  are  established  by  a  rectangle 
with  sides  parallel  to  the  X  and  Y  axes  and  X,  Y  values  of: 

XMAX  =  AMAXI  (XI3,  X23,  X33,  X43) 

XMIN  =  AMINl  (X13,  X23,  X33,  X43) 

YMAX  =  AMAXI  (Y13,  Y23,  Y33,  Yil3) 

YMIN  =  AMINl  (Y13,  Y23,  Y33,  Y43) 

where  all  the  values  in  parentheses  are  subscripted  (IB,  lA).  The  map  point  subscript 
limits  are  then  established  as  IMIN,  IMAX,  3M1N,  JMAX. 

3-3.2  Check  Each  Map  Point  in  the  "Limit  Box"  to  See  if  it  Falls  in  the  Area 

The  approach  used  here  is  the  one  used  by  computer  graphics  routines  to 
determine  whether  a  line  or  point  is  "masked"  by  another  surface.  The  signs  of  the 
cross-products  of  each  of  the  four  sides  of  the  quadrangle  with  the  next  side  srx$2*, 
S2'xS3',  S3'xS4',  and  S4'xSr  are  compared  with  the  signs  of  the  cross-products  of  each 
side  cross  the  vector  between  the  appropriate  quadrangle  vertex  and  the  map  point  in 
question.  This  is  illustrated  in  Figure  3-6,  where  the  sign  of  Sl'xS2*  is  the  same  as  that 
of  ^'xSMPl  (where  map  point  MPl  lies  in  the  quadrangle)  and  is  not  the  same  as 


Sl'xSMP3  (which  lies  outside  the  quadrangle).  However,  the  sign  is  the  same  as  that  of 
SI 'xSMP2  where  MP2  lies  outside  the  quadrangle.  Therefore,  the  signs  of  all  four  sides 
must  be  checked  before  a  map  point  can  be  counted  "in"  the  quadrangle.  However,  the 
map  point  can  be  counted  "out”  of  the  quadrangle  at  the  first  case  of  disagreement  of 
sign.  As  stated  before,  an  additional  check  is  made  to  exclude  points  that  lie  in  the 
"corners"  of  the  quandrangle. 

3-3.3  Increase  Activity  Density  of  Map  Point 

If  the  cross-products  have  agreed  in  sign  in  all  4  cases,  we  can  increase  the 
activity  density  of  the  map  point  1,  3  by  that  of  the  1B,1A  quadrangle: 

DENSMP(I,J)  =  DENSMP(I,a)  +  DENS(1B,IA) 

3-3.#  Check  for  Rough  Agreement  of  Total  Fallout  Activity 

A  check  can  be  made  for  approximate  agreement  between  the  known  total 
fallout  activity  and  the  following  quantity: 

TOT  ACT  =  TOTDEN  *  DEL  *  DEL 

r 

where  TOTDEN  =  DENSMP(I,J) 

This  should  not  necessarily  be  expected  to  yield  close  agreement  since  the  value  of 
DENSMP  (1,3)  is  not  necessarily  uniform  across  its  (1,3)^*^  incremental  area. 

3-4  PLOTTING  THE  MAP  (SUBROUTINE  PLTMAP) 

This  routine  simplifies  the  plotting  procedures  used  in  SEER  because  it  is 
no  longer  necessary  to  plot  shoulders,  edges,  a  separate  GZ  area,  or  do  interpolations. 
The  SEER  map  procedure  has  also  been  changed  to  plot  dots  and  blanks  where  no 
activity  occurs,  as  does  the  DELFIC  map  routine.  The  current  SEER  routine  (CATALG) 
plots  zeros,  and  is  therefore  more  difficult  to  read.  A  print  of  X  and  Y  coordinate 
values  across  the  top  and  down  the  right  side  of  the  map  has  also  been  added.  The 
value  of  TOT  ACT  is  printed  at  the  map  bottom,  but  the  same  caveat  applies  as 
discussed  in  Section  3-3.4. 


3-5  SUMMARY  OF  SUBROUTINES  ADDED  AND  REPLACED 


Subroutines  added  and  replaced  are: 


Subroutines  Added  by  Current  SEER  Subroutines  Replaced  (not  used)  in 

Modified  Approach  Modified  Approach 


RTLFT 

SMEAR 

GRID 

MAPDEN 

PLTMAP 


SEER  3 
SEER  4 
SEER  5 
FODOSE 
CATALG 
DOWN  5 
EDGE5 
GRDZ5 
CIRC5 


CROS5 

DNWD5 

DOSMP 

INTP5 

EANDF5 

EDGY  5 

SHLD5 

UPWD5 


3-6  NEW  PARTICLE  SIZE  GROUP  ACTIVITY  DISTRIBUTION 


For  this  comparison  effort,  SEER  was  modified  to  provide  the  user  with  the 
capability  of  selecting  a  particle  size  group  activity  distribution  which  approximates 
that  currently  used  in  DELFIC.  The  latter  was  obtained  by  modifying  DELFIC  slightly 
(at  GXPSR.6I)  to  print  out  the  percentage  of  activity  of  each  of  its  100  particle  size 
classes.  These  percentages  were  grouped  into  25  classes  to  coincide  as  much  as 
possible  with  the  sizes  of  the  25  SEER  size  classes.  The  user  can  select  this  new 
particle  size  group  activity  distribution  or  the  one  currently  in  SEER  by  the  input  of 
the  value  3  (new)  or  1  (current)  for  the  index  KPGS. 
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SECTION  4 

COMPARISON  RESULTS 


The  comparisons  made  to  judge  the  relative  effectiveness  of  the  two 
versions  of  SEER  were  done  in  two  phases.  The  first  phase  was  done  for  just  one  wind 
and  served  to  establish  the  validity  of  the  new  SEER  methodology.  The  current  SEER 
version  used  did  not  have  the  adjusted  K-factor,  and  thus  was  essentially  not  the 
version  of  SEER  currently  used  in  CIVIC.  The  phase  one  comparison  thus  only 
establishes  the  improvement  of  SEER  over  the  pre-CIVIC  version.  The  phase  two 
comparison  included  the  current  CIVIC  version  of  SEER  with  adjusted  K-factor,  and 
was  done  for  the  multiple  hypothetical  cases  (two  yields,  eleven  winds)  that  were 
treated  in  References  1  and  5. 

4-1  PHASE  ONE  RESULTS  -  COMPARISON  WITH  PRE-CIVIC  SEER 

The  pre-CIVIC  version  of  SEER  (NEW  =  0),  the  modified  SEER  (NEW  =  1), 
and  DELFIC  were  run  for  the  case  of  a  1.0  MT  surface  burst  and  the  winds  observed  at 
0630  25  Dec.  1971  at  Templehof  Airport,  Berlin,  Germany.  Both  SEER  cases  (NEW  = 
(land  NEW  =  1)  were  run  with  the  existing  SEER  particle  size  group  activity 
distribution,  as  well  as  with  an  activity  distribution  which  approximates  that  of 
DELFIC  (see  Section  3-6). 

4-1.1  Comparisons  of  Activity  Areas,  integrated  Activity,  Run  Times 

Table  4-1  summarizes  the  results  obtained  from  running  DELFIC  and  the 
two  versions  of  SEER,  each  with  the  two  particle  size  class  activity  assignments, 
KPGS  =  1  and  KPGS  =  3.  The  following  comments  should  be  noted  in  regard  to  these 
results: 

•  The  areas  enclosed  in  the  activity  contours  were  obtained  in  all  cases 
simply  by  multiplying  the  number  of  map  points  whose  activity  exceeded 
the  particular  level,  by  the  incremental  area  represented  by  each  map 
point  (53.3  km2  for  DELFIC,  103.3  km2  for  SEER  (NEW  =  tf),  and  64.0  km2 
for  SEER  (NEW  =  1). 


Table  4-1.  Summary  of  results  from  DELFIC  and  SEER. 


SEER  (NEW,  KPGS)* 
DELFIC  3  1.  1 

2 

Areas  (in  km  )  Enclosed 
Within  Activity  Contours 
(r/hr  @  1  hr.) 


10 

43,150 

73,680 

74,190 

50,050 

100 

6,133 

9,403 

8,370 

7,552 

1000 

480 

1,343 

826 

704 

3000 

53 

310 

103 

0 

Integ^rated  Activity 


(10^^  r/hr  @  1  hr) 

4.014 

6.796 

5.823 

4.346 

Farthest  Extent  of 

Contour  from  G2  (km) 

10 

709 

705 

713 

725 

100 

201 

255 

275 

215 

1000 

48 

73 

65 

56 

3000 

11 

27 

25 

— 

Run  Time  (sec) 

41.666 

0.654 

0.689 

0.845 

NEW  =  0:  Existing  version  of  SEER  before  modification 

NEW  =  1:  SEER  modified  to  conserve  activity 

KPGS  =  1:  Particle  size  class  activity  assignment  currently  in  SEER 

KPGS  =  3;  Particle  size  class  activity  assignment  from  DELFIC 


•  The  integrated  activity  was  obtained  by  multiplying  the  map  point  activity 
value  by  the  incremental  area  represented  by  each  map  point.  Since  the 
activity  is  not  uniform  from  one  map  point  to  the  next,  the  values  are  not 
expected  to  have  a  high  degree  of  precision. 

•  The  extent  of  the  integrated  activity  summation  was  over  just  the  area  of 
the  fallout  field  that  was  plotted,  and  did  not  include  outlying  areas  of  low 
activity.  Thus,  the  total  value  should  be  expected  to  be  less  than  the 
fallout  portion  of  the  total  activity  within  the  cloud.  Fallout  activity 
values  are  discussed  below  in  Section  4-1.3. 

•  The  furthest  extent  from  the  G2  of  the  4  contour  areas  is  a  program  output 
in  the  DELFIC  program,  and  those  values  are  the  ones  shown.  The  values 
for  SEER  are  measured  from  the  output  maps.  No  values  for  3000  R/HR 
are  given  for  SEER  (NEW  =  1)  because  that  activity  value  was  not  obtained 
for  any  of  the  map  points. 


4-1.2  Contour  Plot  Comparisons 


Figures  4-1  through  4-11  show  comparative  plots  of  activity  contours  for 
10,  100,  1000,  and  3000  r/hr  (at  1  hour  after  detonation).  In  each  figure,  the  DELFIC 
contour  is  compared  to  one  of  the  4  SEER  cases  ((NEW,  KPGS)  =  (0,  1),  (0,  3),  (1,  1), 
(1,  3))  (see  not  at  bottom  of  Table  4-1),  except  that  no  plot  is  shown  for  the  10  r/hr 
case  for  (0,  3)  because  it  does  not  look  significantly  different  than  (0,  1).  The  maps 
plotted  in  this  comparison  for  SEERi  had  map  points  spaced  8  km  apart. 

4-1.3  A  Word  About  Fallout  Activity 

DELFIC  does  not  use  a  K-factor,  or  initialized  value  of  total  activity,  in 
order  to  set  the  level  of  fallout  activity  (p.  26  of  Ref.  4,  p.  8  of  Ref.  8).*  A  fast¬ 
running  code  to  approximate  DELFIC,  such  as  SEER,  must  therefore  start  with  some 
initial  total  activity  level  obtained  from  the  literature.  Table  4-2  shows  values 
provided  in  the  indicated  references.  As  can  be  seen,  there  is  some  disparity  among 
the  values  reported.  Moreover,  there  must  also  be  an  adjustment  for  the  percent  of 
total  activity  that  is  deposited  in  the  so-called  "early  fallout.”  This  early  fallout  is  the 
only  part  that  is  of  strategic  significance  in  the  calculations  to  be  made  with  SEER, 
and  this  fraction  must  therefore  be  multiplied  by  the  total  activity  to  obtain  the  "K- 

factor.”  The  current  SEER  uses  a  K-factor  of  6.0775  x  10^  (r/hr)/(KT/m^)»  which 


♦However,  Reference  10  describes  a  modification  to  DELFIC  to  allow  use  of  a 
K-factor  in  a  conventional  activity  calculation. 
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Figure  4-1.  10  r/hr  ronrour,  DF.LFIC  vs.  old  SEER. 
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Figure  4-2.  10  r/hr  contour,  DELFIC  vs.  modified  SEER. 
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Figure  4-7.  100  r/hr  contour,  DELFIC  vs.  modified  SEER,  new  activity  distribution 
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Table  4-2.  Initial  activity  values  from  the  literature. 


Megacuries 

(R/HR)/(KT/mi2) 

(R/HR)/(KT/m2) 


Ref.  9  Ref.  7  Ref.  6 

450  (p.  5-66)*  530  (p.  453)  - 

-  2900  (p.  390)  3610 

—  7.51  X  10^  9.35  X  10 


%  that  is  included  "roughly  half"  60%  (p.  456) 

in  "early  fallout"  (p.  5-66) 


All  values  are  normalized  to  1  hour  after  detonation  and  are  for 
a  fission  yield  of  1  KT.  The  R/HR  values  are  at  a  height  of  3' 
above  an  infinite  smooth  plane. 


•Reference  9  elaborates  on  this  value  with  the  following  footnote; 

This  value  approximates  the  activity  per  kiloton  in  disintegrations  per 
second  at  one  hour.  Actual  values  may  range  from  about  430  to 
about  460  megacuries,  depending  on  the  fissile  material  and  the 
neutron  spectrum  that  causes  the  fission.  The  user  may  encounter 
other  values  in  various  sources,  e.g.,  550  gamma-megacuries  per 
kiloton  at  1  hour  after  explosion.  This  latter  is  a  fictitious,  but 
useful,  relationship  that  relates  the  fision  product  gamma  source  to 
an  equivalent  monoenergetic  source  with  an  energy  equal  to  the 
average  photon  energy  of  the  fission  products  at  1  hour  after  the 
explosion. 


represents  63%  of  the  Reference  S  activity  value.  DELFIC  output  shows  a  "fraction 
down"  of  0.36S  over  a  wide  range  of  yields. 


9 

For  the  phase  one  comparisons,  the  value  of  6.0775  x  10  was  retained  in 

2 

both  versions  of  SEER.  For  runs  made  to  obtain  the  activity  in  terms  of  curies/m  ,  a 
K-factor  of  31S  megacuries  (60%  of  530)  was  used.  As  discussed  previously,  the 

9 

version  of  SEER  now  being  used  with  CIVIC  hai  a  K  factor  of  only  3  x  10 
2 

(rAir)/(KT/m  ).  This  low  value  was  introduced  in  an  atte  npt  to  produce  lower  results 
and  offset  the  tendency  of  SEER  (NEW  =  0)  to  create  activity.  This  overprediction  of 
activity  has  been  discussed  above,  and  can  be  seen  in  Figures  ^-1,  4-4,  4-5,  4-8,  and 
4-9. 

4-1.4  Discussion  of  Results 

The  phase  one  relative  comparison  of  the  two  versions  of  SEER  (NEW  =  0, 
1)  is  discussed  in  terms  of  the  5  criteria:  (1)  areas  enclosed  in  contours,  (2)  integrated 
activity  over  the  map,  (3)  contour  limit  distances,  (4)  appearance  of  contour  plots,  (5) 
run  times.  The  comparisons  will  be  for  the  SEER  runs  that  used  the  new  (DELFIC) 
particle  group  activity  assignments  (KPGS  =  3).  For  simplicity,  the  old  SEER  will  be 
called  SEERd,  while  the  modified  SEER  will  be  called  SEERl.  (See  Table  4-1  for  the 
values  referred  to  in  this  section.) 

Areas  Enclosed  in  Contours.  SEERl  was  much  closer  to  the  DELFIC  values 
for  the  4  contours  than  was  SEER0.  The  differences  (compared  to  the  DELFIC  areas) 
for  SEERl  were  +14%,  +2%,  -20%,  -100%  (respectively  for  the  10,  100,  1000,  and  3000 
rA>r  contours),  while  for  SEERlI  the  errors  were  +72%,  +36%,  +72%,  +94%. 

Integrated  Activity  Over  the  Map.  The  maps  used  to  plot  the  fallout 
patterns  did  not  cover  the  entire  fallout  in  any  of  the  3  cases.  The  extents  of  the  two 
SEER  maps  were  the  same,  while  the  DELFIC  map  covered  a  larger  area.  However, 
from  the  data  obtained  it  is  clear  that  the  SEERl  integrated  activity  is  closer  to  that 
of  DELFIC  than  is  that  of  SEERd.  The  differences  are  -8%  for  SEERl  and  +45%  for 
SEER0.  Thus,  it  can  be  stated  that,  within  reasonable  limits  of  accuracy,  SEER  1 
appears  to  conserve  activity. 


Contour  Limit  Distances.  SEERl  is  somewhat  better  here,  the  differences 
being  +2%,  +496,  -9%,  -10096;  while  the  differences  for  SEER0  are  +1/2%,  +37%,  +35%, 
+127%,  for  thelO,  100,  1000,  and  3000  r/hr  contours  respectively. 

Appearances  of  the  Contour  Plots.  10  r/hr;  Figure  4-1  clearly  shows  the 
prime  weakness  of  SEER0  in  its  method  of  interpolating,  and  thus  "filling -in"  activity 
where  none  (or  a  lesser  amount)  exists.  SEERl  (Figures  4-2  and  4-3)  does  a  better  job 
of  showing  the  variable  nature  of  the  fallout  field  at  the  outlying  regions,  although  the 
10  r/hr  contour  is  by  no  means  identical  to  that  of  DELFIC.  The  "inner"  SEERl 
contours  shown  in  Figures  4-2  and  4-3  are  actually  "holes"  where  the  dose  rate  is 
slightly  less  than  10  r/hr. 

100  r/hr;  Again  we  see  in  Figures  4-4  and  4-5  the  tendency  of  SEER0  to 
overpredict  activity.  Figures  4-6  and  4-7  show  SEERl  to  resemble  DELFIC  more 
closely. 

1000  r/hr;  The  SEER0  overprediction  is  again  evident  in  Figures  4-8  and 
4-9.  In  Figures  4-10  and  4-11  we  see  that  changing  to  the  DELFIC  particle  size  class 
activity  assignments  dramatically  reduced  the  SEERl  area  coverage  for  the  1000  r/hr 
level  from  a  difference  in  area  from  DELFIC  of  +47%  to  -20%.  The  latter  is  still 
closer  to  DELFIC  than  the  +72%  of  SEER0. 


3000  r/hr;  SEERl  did  not  produce  a  map  point  at  this  level  of  dose  rate. 
While  SEER0  did  produce  1  map  point  at  3000  r/hr  and  2  others  above  2500,  the 
location  does  not  coincide  with  that  of  DELFIC. 

Run  Times.  The  run  times  for  SEER0  and  SEERl  are  comparable,  with  the 
latter  being  19%  longer  for  KPGS  =  3.  However,  run  times  were  seen  to  vary  from  one 
run  to  another,  and  no  attempt  was  made  to  streamline  the  final  version  of  SEERl  or 
make  a  methodical  comparison  over  several  runs  of  the  "final"  version  of  each 
(SEER0  and  SEERl).  Most  of  the  unnecessary  diagnostic  write  statements  were 
removed,  however,  and  the  results  shown  in  Table  4-1  are  sufficiently  representative 
to  state  that  SEERl  increases  the  run  time  over  SEER0  by  no  more  than  20%  or  25% 
on  the  average. 


♦-2  PHASE  TWO  RESULTS  -  COMPARISONS  WITH  CIVIC  VERSION  OF  SEER 

Following  the  initial  comparisons  of  Phase  1,  a  more  comprehensive  set  of 
comparisons  was  made  which  included  the  current  CIVIC  version  of  SEER  (with  the 
adjusted  K-factor).  A  set  of  previously  run,  hypothetical  shot,  cases  was  available  for 
use  in  these  comparisons  (References  1  and  5).  These  cases  covered  11  different  winds 
and  2  yields  (50  KT  and  1.0  MT).  The  results  that  will  be  discussed  here  are  for  the  1.0 
MT  hypothetical  shot  yield  only,  because  the  50  KT  results  were  essentially  the  same 
and  lead  to  the  same  conclusions.  Because  of  the  data  available  from  References  1 
and  5,  the  results  were  compared  with  regard  to  areas  enclosed  in  the  contours  (100, 
300,  1000,  and  3000  r/hr)  and  contour  limit  distances  only.  Contour  areas  only  are 
shown  here  because  the  results  and  conclusions  are  the  same  for  contour  limit 
distances. 


The  plots  that  are  shown  in  Figures  4-12  through  4-16  show  the  contour 
area  results  in  square  kilometers  for  all  eleven  winds  and  for  the  four  dose  rate 
contours.  SEER  contour  areas  are  plotted  vs.  those  of  DELFIC  (SEER  along  the 
ordinate  and  DELFIC  along  the  abscissa).  Areas  that  are  equal  (if  any)  would  be 
plotted  along  the  45°  line  shown.  All  plots  shown  are  for  1.0  MT  yield. 

From  comparison  of  Figures  4-12  and  4-13  it  can  be  seen  that  the  existing 
version  of  SEER  comes  into  much  better  agreement  with  DELFIC  when  the  adjusted 
K-factor  is  applied.  When  the  particle  size  group  activity  distribution  is  changed  to 
approximate  that  of  DELFIC  (Figure  4-14),  the  results  depend  upon  the  dose-rate 
contour.  The  SEER-DELFIC  correlation  is  slightly  worse  at  the  100  and  300  r/hr 
contours,  about  the  same  for  1000  r/hr,  and  somewhat  better  for  the  3000  r/hr 
contour. 


The  modified  SEER  (Figure  4-15)  provides  a  dramatic  improvement  over 
the  old  (pre-CIVIC,  no  adjustment  to  K-factor)  SEER,  but  is  superior  to  SEER/CIVIC 
(adjusted  K-factor)  only  for  the  100  and  3000  r/hr  contours  (this  superiority  is 
eliminated  at  3000  r/hr  by  use  of  the  DELFIC  particle  size  group  activity  distribution 
for  SEER/CIVIC).  Application  of  the  DELFIC  particle  size  group  activity  distribution 
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Figure  4-15.  Contour  areas,  DELFIC  vs.  modified  SEER. 
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to  the  modified  SEER  (Figure  4-16)  worsens  the  approximation  to  DELFIC  at  the  1000 
and  3000  r/hr  contours,  and  therefore  accentuates  the  fact  that  SEER  lacks  a  model 
for  the  cloud  stem. 


SECTION  5 


CONCLUSIONS 


It  is  apparent  that  the  modification  to  SEER  reported  here  has  corrected 
one  deficiency  of  SEER  (non-conservation  of  activity),  but  has  served  also  to 
accentuate  the  other  (lack  of  a  cloud  stem).  The  modified  SEER  has  a  more  accurate 
physical  basis  than  the  current  SEER,  and  therefore  is  a  more  faithful  representation 
of  activity  deposition,  given  the  debris  cloud  modeled  by  the  code.  Additional  work 
would  have  to  be  done  to  overcome  the  no-stem  deficiency.  The  additional  run  time  of 
the  modified  SEER  is  not  considered  critical,  but  great  care  would  have  to  be 
exercised  in  addressing  the  stem  deficiency  to  avoid  significant  increases  to  run  time. 
The  overall  conclusion  of  this  effort  (and  considering  the  results  reported  in  Reference 
1)  is  that,  until  further  work  can  be  done  on  the  modification,  the  current  CIVIC 
version  of  SEER  with  adjusted  K-factor  is  the  best  available  fast-running 
approximation  to  DELFIC. 


*  *  1  ' 
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